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ABSTRACT 

A  MACSYMA  program  is  presented  which  determines  whether  a  given  single  nonlinear  ODE  or  PDE 
with  (real)  polynomial  terms  fulfills  the  necessary  conditions  for  having  the  Painleve  property.  Together 
with  some  mathematical  background,  we  give  a  synopsis  of  the  algorithm  for  the  program,  its  scope  and 
limitations.  Various  examples  of  typical  output  of  the  program  are  provided. 
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3  Introduction:  Mathematical  Background 

Roughly  speaking  dynamical  systems  may  be  divided  into  two  classes.  In  the 
first  class  one  finds  those  systems  which  exhibit  chaotic  behavior,  i.e.  the  solu¬ 
tions  of  such  systems  depend  sensitively  on  the  initial  data.  A  consequence  of 
this  chaotic  behaviour  is  that  such  systems  are  usually  not  explicitely  solvable 
in  terms  of  “elementary  functions.”  The  other  class  of  systems  contains  those 
equations  which  are  algebraically  completely  integrable.  Obviously,  one  wants  to 
characterize  systems  which  do  not  lead  to  chaos. 

At  the  turn  of  the  century,  Pain' we  and  his  collaborators  were  able  to  iden¬ 
tify  all  second  order  ODEs  of  the  form  fzz  =  K(z,  /,  fz),  which  are  globally 
integrable  in  terms  of  elementary  functions  by  quadratures  or  by  linearization. 
The  restrictions  on  the  function  K,  which  is  rational  in  fXl  algebraic  in  /,  and  an¬ 
alytic  in  x,  arise  from  careful  singularity  analysis.  Indeed,  integrability  requires 
that  the  only  movable  singularities  in  the  solution  f{x)  are  poles.  Singulari¬ 
ties  are  movable  if  their  location  depends  on  the  initial  conditions.  Hence,  the 
critical  points  (including  logarithmic  branch  points  and  essential  singularities) 
ought  to  be  fixed  to  have  integrability  and  consequently  to  assure  predictable 
behavior  of  solutions  for  all  initial  conditions. 

Definition:  A  simple  equation  or  system  is  said  to  have  the  Patnleve  property 
(PP)  if  its  solution  in  the  complex  plane  has  no  worse  singularities  than  movable 
poles. 

In  a  broader,  more  modern  context,  including  PDEs,  integrability  became 
associated  with  the  existence  of  a  Lax  representation  which  allows  linearization 
of  the  given  equation(s).  The  solution  may  then  be  constructed  by  e.g.  the 
Inverse  Scattering  Transform. 

A  restricted  class  of  completely  integrable  PDEs,  admits  solitary  traveling 
wave  solutions  (called  soliions  if  they  conserve  their  identity  upon  collision). 
Completely  integrable  systems  of  ODEs  admit  a  finite  dimensional  Hamiltonian 
formulation  and  have  a  finite  number  of  first  integrals.  Analogously,  completely 
integrable  PDEs  possess  an  infinite  number  of  conserved  quantities,  infinitely 
many  symmetries,  nontrivial  prolongation  structures  and  Kac-Moody  algebras. 

1  This  work  is  partially  supported  by  AFOSR  under  Grant  85-NM-0263.  We  thank 
Jonathan  Len  (Symbolics,  Inc.)  for  writing  the  POWERS  subroutine.  The  original  co¬ 
developer  of  the  program,  Dr.  Eric  Van  den  Bulck,  presently  at  the  University  of  Leuven 
(Belgium),  is  gratefully  acknowledged  for  endless  nights  of  debugging.  Help  from  Dr.  D. 
Rand  (University  of  Montreal)  in  further  debugging  of  the  program  was  very  much  appreci¬ 
ated. 
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In  a  first  attempt  to  identify  integrable  PDEs,  Ablowitz  et  al.  [for  refer¬ 
ences:  consult  the  books  and  review  papers  cited  at  the  end]  conjectured  that 
every  ODE,  obtained  by  an  exact  reduction  of  a  PDE  solvable  by  1ST,  possesses 
the  Painleve  property.  This  conjecture  was  based  on  the  fact  that  well-known 
integrable  PDEs,  such  as  the  Korteweg-de  Vries  (KdV)  equation,  can  be  trans¬ 
formed  (by  similarity  transformations)  intoODEs  of  Painleve  type,  in  particular 
into  the  so-called  six  Painleve  transcendents.  This  is  yet  another  indication  that 
integrability  is  closely  related  to  the  absence  of  movable  critical  points. 

2  Algorithm 

An  algorithm,  due  to  Weiss  ei  al  enables  us  to  verify  if  the  given  ODE  or  PDE 
satisfies  the  necessary  criteria  to  have  the  Painleve'  property. 

For  the  PDE  case:  The  solution  /,  say  in  only  two  independent  variables 
(<,  x),  expressed  as  a  Laurent  series, 

OO 

f  =  gaY^Uksk  u) 

k  =  0 

should  only  have  movable  poles.  In  (1),  u0(t,x)  ^  0,  a  a  negative  integer, 
and  Ufc(t.x)  are  analytic  functions  in  a  neighborhood  of  the  singular;  non- 
characteristic  manifold  g{t,x)  =  0,  with  pr(t,x)  ^  0. 

For  the  ODE  case2:  x  will  be  replaced  by  g  +  x0  in  (1);  xo  being  the  initial 
value  for  x. 

The  Painleve'  test  is  carried  out  in  three  steps: 

•  Step  1:  Determine  the  negative  integer  a  and  «o  from  the  leading  order 
“ansatz”  /  oc  Uoga ,  by  balancing  the  minimal  power  terms  after  substitu¬ 
tion  of  this  first  term  of  (1)  into  the  given  equation 

•  Step  2:  Calculate  the  non-negative  1.  r  powers  r,  called  the  reso¬ 
nances,  at  which  arbitrary  functions  uT  ■  n(  .r  ?he  expansion.  This  is  done 
by  requiring  that  ur  is  arbitrary  after  su  ution  of  /  oc  u0g°  +  urg°+r 
into  the  equation,  only  retaining  its  most  singular  terms 

•  Step  3:  Verify  that  the  correct  number  of  arbitrary  functions  ur  indeed 
exists  by  substituting  a  truncated  expansion  of  the  form  (1),  in  which  one 
sums  over  I:  =  1,2,  ...,rmax,  where  rmax  represents  the  largest  resonance, 
into  the  given  equation. 

At  non-resonance  levels,  determine  ut  unambiguously. 

At  resonance  levels,  ur  should  be  arbitrary  due  to  a  vanishing  coefficient  of 
gr+minpowg  One  check  whether  or  not  this  compatibility  condition 

JFor  this  case  a  MACSYMA  package  already  existed  [Ran] 
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is  satisfied.  Here  minpowg  denotes  the  (negative)  power  in  g  of  the  most 
singular  terms  in  the  equation. 

An  equation  or  system  for  which  the  above  steps  can  be  carried  out  consis¬ 
tently,  and  for  which  the  compatibility  conditions  at  all  resonances  are  satisfied, 
is  said  to  have  the  Painleve  Property  and  is  conjectured  to  be  integrable. 

The  reader  should  be  warned  that  the  above  algorithm  does  not  detect  the 
existence  of  essential  singularities.  In  other  words,  for  an  equation  to  be  inte¬ 
grable  it  is  necessary  but  not  yet  sufficient  that  it  passes  the  Painleve  test. 
There  are  integrable  equations  (see  Ex.  2)  that  only  have  the  Painleve'  property 
after  a  suitable  chance  of  variables. 

Refinements  of  the  Painleve  property  have  been  established,  allowing  for 
rational  power  expansions  of  /,  and  hence  including  certain  algebraic  branch 
points  in  addition  to  movable  poles. 

3  Scope  and  Limitations  of  the  Program 

3.1  Scope 

•  the  program  works  for  a  single  ODE  or  PDE 

•  the  degree  of  nonlinearity  in  all  the  variables  is  unlimited 

•  the  number  of  parameters  in  the  equation  is  unlimited 

•  the  number  of  independent  variables  is  also  unlimited3 

•  ODEs  and  PDEs  may  have  explicitly  given  space  (and/or  time)  dependent 
coefficients  of  integer  degree  (see  Ex.  4) 

•  PDEs  may  have  arbitrary  time  and  space  dependent  coefficients  (see  Ex. 

4) 

•  coefficients  may  be  complex,  although  the  usefulness  of  the  Painleve  test 
is  then  debatable 

•  a  selected  positive  or  negative  rational  value  of  a,  or  a  =  0  can  be  supplied 
by  the  user 

•  the  time  consuming  calculation  of  the  coefficients  uu  and  the  verification 
of  the  compatibility  conditions  is  made  optional 

•  it  is  possible  to  substitute  an  expansion  of  the  form  (1)  with  a  selected 
number  of  terms,  e.g.  to  carry  on  with  the  calculations  beyond  rmax 

3The  released  version  of  the  program  works  for  at  most  four  variables  ((, r,  y,  r)  but  this 
can  be  extended  trivially 
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•  the  output  provides  vital  information,  including  error  messages  and  warn¬ 
ings,  to  remedy  possible  problems 

3.2  Limitations 

•  systems  of  equations  are  presently  excluded 

•  the  test  is  restricted  to  the  traditional  Painleve  test  based  on  the  expansion 
(1),  with  at  least  rational  a,  hence  general  fractional  expansions  in  rj  arc 
at  present  not  possible 

•  transcendental  terms  in  the  equation  are  not  allowed,  but  they  can  often 
be  removed  by  a  suitable  (exponential)  transformation  of  the  dependent 
variable  (see  Ex.  2) 

•  arbitrary  parameters  in  the  powers  of  /  and  its  derivatives  are  not  allowed 

•  neither  are  arbitrary  (unspecified)  functions  of  /  and  its  derivatives 

•  selective  substitution  of  certain  u*  is  presently  not  possible,  i.e.  u0,  ult 
etc.  are  explicitly  determined  whenever  possible,  and  their  expressions  are 
used  in  the  calculation  of  the  next  ut 

•  apart  from  trivial  cases,  nonlinear  equations  for  uq  are  not  solved,  if  they 
occur  the  program  carries  on  with  the  undetermined  coefficient  u0  (see 
Ex.  5) 

•  the  program  only  checks  whether  or  not  the  compatibility  condition  is 
satisfied.  It  does  not  solve  for  arbitrary  parameters  (or  functions)  or  for 
«o  and  its  derivatives,  should  these  occur  (see  Exs.  3,  4  and  5) 

•  intermediate  output  is  only  possible  by  including  extra  print  statements 
in  the  source  code  of  the  program 

•  the  expressions  occuring  in  the  output  on  the  screen  are  not  accessible  for 
further  interactive  calculations  without  modification  of  the  program 
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4  Spin-offs  of  the  Painleve  Analysis 

The  Painleve  test,  apart  from  it  usefulness  in  testing  the  integrability  of  Hamil¬ 
tonian  systems,  evolution  and  wave  equations,  has  rather  interesting  connections 
with  standard  techniques  in  the  study  of  dynamical  systems  and  soliton  theory: 

•  Truncation  of  the  Laurent  series  (1)  at  the  constant  level  term  leads  to 
auto-Backlund  transformations 

•  The  resulting  Painleve-Backlund  equations,  obtained  by  substitution  of 
the  truncated  expansion  and  equating  to  zero  powers  of  g,  can  readily  be 
linearized  to  derive  Lax  pairs  for  various  systems  of  ODEs  and  PDEs 

•  As  a  consequer.-e  for  ODEs,  it  is  possible  to  construct  algebraic  curves 
and  explicitly  integrate  the  equations  of  motion 

•  For  PDEs,  Painleve  analysis  provides  insight  in  the  construction  of  single 
and  multi-soliton  solutions  and  also  rational  solutions,  via  direct  methods 
(Hirota’s  formalism  and  its  clones)  (see  Exs.  3  and  5) 

•  The  Painleve'  property  serves  as  in  tool  in  identifying  the  infinite  dimen¬ 
sional  symmetry  algebras  for  PDEs,  which  in  turn  have  the  structure  of 
subalgebras  of  Kac-Moody-Virasoro  algebras 

5  Using  the  Program 

The  program  carries  out  the  Painleve  test  in  batch  mode  without  interaction 
by  the  user.  The  user  only  has  to  type  in  the  LHS  of  the  equation  (denoted  by 
eq)  and  possibly  select  some  options. 

•  For  ODEs:  The  use  of  dependent  variable  /  and  independent  variable  x 
is  mandatory.  A  typical  term  in  the  ODE  reads  fx[.]( x),  where  within 
the  brackets  the  order  of  derivation  is  inserted.  The  function  without 
derivatives  may  be  denoted  by  /  itself.  The  symbol  eq  denotes  the  L11S 
of  the  equation 

•  Ex.:  To  test  the  Fisher  ODE,  frI  +  afr  —  f2  +  /  =  0,  one  would  enter 

eq  :  /x[2](r)  +  o  *  /[x](x)  +  /  *  *  2  -  /; 

The  program  will  then  treat  a  as  an  arbitrary  parameter. 

•  For  PDEs:  Analoguously,  a  typical  term  reads  ftxyz[k,  l,  m,  n](t,  x,  y,  z), 
where  the  integers  k,l,m ,  and  n  are  the  orders  of  derivation  with  respect 
to  the  variables  <,x,y,  and  z 
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•  Ex.:  To  test  the  KdV  equation,  }t  +  affx  +  fxtI  =  0,  one  enters 
eq  :  ftx[l,  0](t,  x)  +  a*  f  *  ftx[  0,  l](f,  ar)  +  ftx[  0,  3](t,i); 

Again,  the  program  will  treat  a  as  an  arbitrary  parameter. 

The  program  is  available  (on  IBM  floppy  or  Macintosh  disk)  for  two  slightly 
different  version  of  MACSYMA:  Franz  Lisp  MACSYMA  309  and  Common  Lisp 
MACSYMA  412.  Tree  copies  may  be  obtained  at  the  above  address.  To  get 
the  program  by  email,  contact:  HEREMAN@EUNICE.MRC.WISC.EDU  (arpa 
net). 

The  program  has  been  successfully  tested,  with  MACSYMA  versions  309.3 
and  412,  on  a  VAX  2000  and  a  VAX  11/780,  for  more  than  100  ODEs  and 
PDEs. 

©Willy  Hereman  and  Eric  Van  den  Bulck:  No  part  of  the  program  PAINLEVE 
TEST  may  be  reproduced  or  sold  without  written  consent  of  the  authors. 


6  Examples 

The  selected  examples  reflect  some  of  the  capabilities  of  the  program.  In  the  five 
examples,  a,  b  and  c  are  arbitrary  parameters,  and  a(t)  is  an  arbitrary  function. 
Many  more  examples,  tested  with  this  program,  may  be  found  in  the  cited  books 
[Lea,  Ste]. 

Example  1:  The  Korteweg-de  Vries  Equation 

For  the  ubiquitous  KdV  equation, 

ft  +  ffx  +  A/rrr  =  0,  (2) 

the  program  provides  the  following  output: 


PAINLEVE  ANALYSIS  OF  EQUATION,  bfxxz  +  ffz  +  ft  =  0 

SUBSTITUTE  u0  ga!fa  FOR  f  IN  ORIGINAL  EQUATION. 
MINIMUM  POWERS  OF  g  ARE  [2  alfa  -  1  ,a//a  -  3] 

*  COEFFICIENT  OF  g2  al'a~l  IS  u0 2  alfa  gx 

*  COEFFICIENT  OF  gal^~3  IS  u0  ( alfa  -  2)  ( alfa  -  1)  alfa  b  (<?*)3 


FOR  EXPONENTS  (2  alfa  -  1)  AND  ( alfa  -  3)  OF  g,  DO 
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WITH  alfa  =  -2,  POWER  OF  g  is  -5  — .  SOLVE  FOR  u0 
TERM  -2  u0  gx  (12  b  ( gx )2  +  u0)  jr  IS  DOMINANT 
IN  EQUATION. 

WITH  u0  =  -12  b  ( gz )2  — .  FIND  RESONANCES 
SUBSTITUTE  u0  galla  +  ur  gr+al}a  FOR  f  IN  EQUATION 

TERM  b  (gz)3  (r  -  6)  (r  -  4)  (r  +  1)  ur  gr~ 5  IS  DOMINANT 
IN  EQUATION. 

THE  2  NON-NEGATIVE  INTEGRAL  ROOTS  ARE  [r  =  4,  r  =  6] 
WITH  MAXIMUM  RESONANCE  =  6  — ♦  CHECK  RESONANCES. 
SUBSTITUTE  POWER  SERIES  £®=o  ut  FOR  f  IN  EQUATION. 
WITH  u0  =  -12  6  (Sx)2 

*  COEFFICIENT  OF  T-  IS  6b(gz):!((-12b(gz)2)x  -  36bgzgxz  +  5 ul3z) 

«i  —  126<yxx 

*  COEFFICIENT  OF  ±  IS  24bgz(4bgzgzzz -3b(gzz)2 +  u2(gz)2  +  gtgz) 

_  4bgzgzzz  -  3b(g  rr)2 

u2  —  /  v  2 

(?r)‘ 

*  COEFFICIENT  OF  -L  IS 

9 

—  \2b(b(gz)2gxzzz  -  4bgxgxzgzzz  +  36(Sxx)3  -  gtgzgzz 

-«3(ffr)4  +  gtx{9z)i)/9x 

_  b(gz)  gzxzx  —  4bgzgzzgXXI  Hh  36(gxx)  —  gigxgZr  4-  gn(gz') 

*  COEFFICIENT  OF  1  IS  0 

3 

u4  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION  IS  SATISFIED  ! 

*  COEFFICIENT  OF  !  IS 

—  —  9b3(gz)  gxx9xxzzz  —  17fc^(jr)^3rrr?rrrr 

+4862(5x)2(</xx)2jxiix  -  26<7i(9r)3<7xxxx  +  7062(yx)2^xx(<7xxx)2 

—  17462yx(jxx)3yxxr  +  l7bgt(gz)2gxZgxxx  -  86<?<r(<7x)3</rrr 
+8l62(<7xt)5  -  2l60t<7x(<7tx)3  +  2l6p,x(<7x)2(5xx)2 

+61146(3!^  gxx  —  965ixx(9i)  Str  +  (st)  (Sx)  3xx  4"  6u56(<;x)® 
+6(u4)t6(5x)7  +  gtt{gx)*  +  26j?,xxi(<7x)4  -  2 0<0<r(ffx)3)/(ffx)5 
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t»5  = 


-(fc2(0r)40rrrrn  “  962(pr )3<?rrSxrrxr  ~  l~b2(gz)3  OiTrOzrrr 

+486 2(,gzf{gzz)2gzzzz  -  2bgt(gz)3gzzzz  +  10b2(gz)2 gzz(gzzz)2 
—  1 1 46* gz{9zz)  9zzz  +  1  •'  bgt{gt)~ gzzgzzz  ~  Sbgtz{gz )■  gzzz 
+8162(pn)5  -  2lbgtgz(gzx)3  +  21bgtx{gz)2 {gzz)2  +  QuAb(gz )6gzz 
-9bgtzz(gz)3gzz  +  {g,)2{gz)2 gzz  +  6(u4)r6(yr)7  +  gtt(gz j1 
+2bgtzzz(gzy  -  2gtgtz{gz)3)/(&b{gz)6) 

*  COEFFICIENT  OF  g  IS  0 
u6  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION  IS  SATISFIED  ! 


Example  2:  The  sine-Gordon  Equation 

Tlie  transcendental  term  in  the  sine-Gordon  equation,  in  light  cone  coordinates. 

u(r  —  sin(u)  =  0,  (3) 

can  be  removed  by  the  simple  substitution  /  =  exp(iu)  to  obtain  an  equivalent 
equation  with  polynomial  terms: 

—2ftfz  +  2f  ftz  —  f3  +  f  =  0.  (-1) 


PAINLEVE  ANALYSIS  OF  EQUATION,  -2 f,fz  +  2 ffiz  -  f3  +  f  =  0 

SUBSTITUTE  u0  ga,>a  FOR  f  IN  ORIGINAL  EQUATION 
MINIMUM  POWERS  OF  g  ARE  [2a//a  -  2,  3  alfa,  alfa ] 

*  COEFFICIENT  OF  g2  a,la~2  IS  -2  u02  alfa  gt  gz 

*  COEFFICIENT  OF  g3a,Ia  IS  -u03 

*  COEFFICIENT  OF  ga,J°  IS  u0 


FOR  EXPONENTS  2  at  fa  -  2  AND  3  alfa  OF  g,  DO 

WITH  alfa  =  -2,  POWER  OF  g  is  6  — >  SOLVE  FOR  u0 
TERM  uo2  (4g,gz  -  u0)±  IS  DOMINANT  IN  EQUATION. 
WITH  u0  =  4gtgz  — *  FIND  RESONANCES 
SUBSTITUTE  u0  9a,,a  +  ur  gr+a'la  FOR  f  IN  EQUATION 


TERM  8(<7«)2($r)2(r  -  2)(r  +  1)  ur  gr~6  IS  DOMINANT 
IN  EQUATION. 

THE  ONLY  NON-NEGATIVE  INTEGRAL  ROOT  IS  [r  =  2] 

WITH  MAXIMUM  RESONANCE  =  2  — *  CHECK  RESONANCES. 
SUBSTITUTE  POWER  SERIES  £LoS*“2«*  FOR  f  IN  EQUATION. 
WITH  u0  =  4 gtgz 

*  COEFFICIENT  OF  ±  IS  -16(5l)2(4s<r  +  ui)(jr)2 
«i  =  -4  gtx 

*  COEFFICIENT  OF  ^  IS  0 
u2  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION  IS  SATISFIED  ! 


FOR  EXPONENTS  (2  alfa  -  2)  AND  {alfa)  OF  g ,  alfa  =  2  IS 
NON-NEGATIVE. 


FOR  EXPONENTS  (3  alfa)  AND  {alfa)  OF  g ,  alfa  =  0  IS  NON-NEGATIVE. 


Example  3:  The  Fisher  Equation 

^From  rigorous  analysis  it  follows  that  if  the  initial  datum  is  given  by  u(0,  z)  =  1 
(z  <  0),  u(0,  z)  =  0  (z  >  0),  then  the  solution  of  the  Fisher  equation, 

ti(  -  azz  +  u2  -  u  =  0,  (5) 

will  converge  to  a  travelling  wave  of  speed  c  =  2.  Furthermore,  for  every  speed 
c  >  2  there  is  a  travelling  wave  with  u(<,-oo)  =  l,u(t,oo)  =  0.  In  1979,  an 
exact  closed  form  solution  of  (5)  was  constructed:  u{t,x)  =  U (z  —  zo  — 
where 


with  z0  any  constant.  The  Painleve  analysis  for  (5),  put  into  a  travelling  frame  of 
reference,  exactly  determines  this  particular  wave  speed  c  =  which,  indeed, 
is  larger  than  2. 


PAINLEVE  ANALYSIS  OF  EQUATION,  /„  +  cfz  -  /2  +  /  =  0 


SUBSTITUTE  u0  ga,fa  FOR  f  IN  ORIGINAL  EQUATION. 
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MINIMUM  POWERS  OF  g  ARE  [2  alfa,  alfa  -  2] 

*  COEFFICIENT  OF  g'2a,la  IS  -u02 

*  COEFFICIENT  OF  g“'l°-*  IS  u0  {alfa  -  1)  alfa 


FOR  EXPONENTS  (2  alfa)  AND  {alfa  -  2)  OF  g,  DO 

WITH  alfa  -  -2,  POWER  OF  g  is  -4  — .  SOLVE  FOR  u0 
TERM  ~(u0  -  6)  u0  -r  IS  DOMINANT  IN  EQUATION. 

WITH  ti0  =  6  — .  FIND  RESONANCES 
SUBSTITUTE  u0  ga,Ja  +  ur  gr+a'J“  FOR  f  IN  EQUATION 
TERM  (r  -  6)(r  +  1)  ur  gr~ 4  IS  DOMINANT  IN  EQUATION. 

THE  ONLY  NON-NEGATIVE  INTEGRAL  ROOT  IS  [r  =  6] 

WITH  MAXIMUM  RESONANCE  =  6  — »  CHECK  RESONANCES. 
SUBSTITUTE  POWER  SERIES  Zt=0gk~2nk  FOR  f  IN  EQUATION. 
WITH  u0  =  6 


*  COEFFICIENT  OF  ±  IS  -2  (6c  +  5m) 


« l  =  -- 


6c 


*  COEFFICIENT  OF  \  IS  - 


6(c2  +  50u2  -  25) 
25 


ti2  ~  -- 


(c  -  5)(c  +  5) 
~50 


*  COEFFICIENT  OF  1  IS 


6(c3  +  25OU3) 
125 


« 3  =  - 


250 


*  COEFFICIENT  OF  1  IS  - 
7c4  -  125 


u  4  =  -  ■ 


5000 


*  COEFFICIENT  OF  5  IS  - 
I3"5) 


c(79c4 
u5  = - 


7c4  +  5000u<  -  125 
500 


79c5  -  1375c +  75000ii5 
12500 


75000 
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*  COEFFICIENT  OF  g 2  IS  - 


c2(6c2  —  25)(6c2  +  25) 
6250 


=  0 


u6  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION:  —  =  0, 


*  *  *  CONDITION  IS  NOT  SATISFIED  .*  *  * 

***  CHECK  FOR  FREE  PARAMETERS  OR  PRESENCE  OF  u0  *** 


Example  4:  The  cylindrical  KDV  Equation 
The  cylindrical  Korteweg-de  Vries  equation, 

~  +  fzxzt  +  6//rI  +  6(/r)2  -f  ftx  =  0,  (7) 

has  the  Painleve  property.  One  easily  determines  the  coefficient  ^  in  (7),  by 
analyzing  a  cylindrical  KdV  equation  with  arbitrary  coefficient  a(t)  of  fT.  Inte¬ 
gration  of  the  compatibility  condition  a(£),  +  2a(£)2  =  0,  gives  a(t)  = 


PAINLEVE  ANALYSIS  OF  EQUATION, 

«(0/x  4"  fzzzx  +  6ffxx+6{fx)7  +  ftx  =  o 


SUBSTITUTE  u0  ga,fa  FOR  f  IN  ORIGINAL  EQUATION. 

MINIMUM  POWERS  OF  g  ARE  [2  alfa  -  2,  alfa  -  4] 

*  COEFFICIENT  OF  g 2  a^a~7  IS  6u02  alfa  (2 alfa  -  1)  (<?r)2 

*  COEFFICIENT  OF  galf°~A  IS  u0  ( alfa-3)(alfa-2)(alfa-l )  a//a  (?x)4 


FOR  EXPONENTS  (2  a//a  -  2)  AND  (a//a  -  4)  OF  g,  DO 
WITH  a//a  =  -2,  POWER  OF  g  is  -6  — -  SOLVE  FOR  u0 
TERM  60  u0  {gzy  (2 (gx)2  +  u0)  ^  IS  DOMINANT 
IN  EQUATION. 

WITH  uo  =  -2(<?r)2  — *  FIND  RESONANCES 
SUBSTITUTE  u0  gal>a  4-  ur  gr+alfa  FOR  f  IN  EQUATION 
TERM  (?r)V  -  6)(r  -  5)(r  -  4)(r  +  l)ur<f-6  IS  DOMINANT 
IN  EQUATION. 
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THE  3  NON-NEGATIVE  INTEGRAL  ROOTS  ARE 
[r  =  4,  r  =  5,  r  =  6] 

WITH  MAXIMUM  RESONANCE  =  6  — -  CHECK  RESONANCES. 
SUBSTITUTE  POWER  SERIES  jfts.0gk~3uk  FOR  f  IN  EQUATION. 
WITH  tt0  =  -2 (gx)2 

*  COEFFICIENT  OF  T.  IS  120(ffr)‘,(25rr  -  ut) 

ui  "  2 gzz 

*  COEFFICIENT  OF  \  IS 

a ' 

~12(9x)~(4gx9zxx  ~  3(yri)"  +  6u2(gx)~  +  9tgz) 

a9x  9xxz  -  3(gtr)2  +9tgt 
2  6  (gzy 

*  COEFFICIENT  OF  -L  IS 

a3 

4((<Ir)  a(0  +  (<7r)2<7rrrr  —  4 9z9zz9xzz 

+3(ffn)2  -  9tgz9xx  -  6 u3(gzy  +  ffn(ffx)2) 


u3  —  ((?x)  a(0  +  (?r)2<?rrrr  ^9z9xz9xzz  +  3(^xr)"  —  9i9z9zz 

+9ix(9x)2)/{G{9x)4) 

*  COEFFICIENT  OF  -L  IS  0 
a 2 

«4  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION  IS  SATISFIED  ! 

*  COEFFICIENT  OF  ±  IS  0 

a 

«5  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION  IS  SATISFIED  ! 

*  COEFFICIENT  OF  1  IS 

6 

vq  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION:  -  0> 

*  *  *  CONDITION  IS  NOT  SATISFIED  .*  *  * 

**«  CHECK  FOR  FREE  PARAMETERS  OR  PRESENCE  OF  u0  »** 
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Example  5:  The  Fitzhugh-Nagumo  Equation 

We  recently  discovered  that  the  Fitzhugh-Nagumo  equation, 

-  ti„  -  u(l  -  u)(u  -  a)  =  0,  (8) 

has  a  closed  form  travelling  wave  solution,  u(t,  r)  =  U(x  —  ct),  where 


u(  0  = 


(9) 


and  c  =  Carrying  out  the  Painleve  test  for  (8),  in  a  travelling  frame,  leads 

to  a  compatibility  condition  which  for  uq  =  \/2  factor  into 


So,  we  see  that  one  of  the  roots  corresponds  with  the  wave  speed  in  (9). 

PAINLEVE  ANALYSIS  OF  EQUATION,  /rr  +  c/t-/3  +  (a+l)/2-a/  =  0 


SUBSTITUTE  u0  ga,fa  FOR  f  IN  ORIGINAL  EQUATION. 
MINIMUM  POWERS  OF  g  ARE  [3a//a,  2 alfa,  alfa  -  2] 

*  COEFFICIENT  OF  g3al'a  IS  -u03 

*  COEFFICIENT  OF  g2al!a  IS  u02  (o  +  1) 

*  COEFFICIENT  OF  gal>a~ 2  IS  u0  ( alfa  -  1)  alfa 


FOR  EXPONENTS  (3 alfa)  AND  (2 alfa)  OF  g,  alfa  =  0  IS  NON-NEGATIVE. 


FOR  EXPONENTS  (3 alfa)  AND  ( alfa  -  2)  OF  g,  DO 

WITH  alfa  =  -1,  POWER  OF  g  is  -3  — -  SOLVE  FOR  u0 
TERM  -u0(u02  -  2)js  IS  DOMINANT  IN  EQUATION. 

WITH  u02  =  2  — >  FIND  RESONANCES 
SUBSTITUTE  u0  ga,*a  +  ur  gr+alfa  FOR  f  IN  EQUATION 
TERM  (r  -  4)(r  +  1)  ur  gr~3  IS  DOMINANT 
IN  EQUATION. 

THE  ONLY  NON-NEGATIVE  INTEGRAL  ROOT  IS  [r  =  4] 
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WITH  MAXIMUM  RESONANCE  =  4  — >  CHECK  RESONANCES. 
SUBSTITUTE  POWER  SERIES  FOR  f  IN  EQUATION. 

WITH  u02  =  2 

*  COEFFICIENT  OF  ±  IS  -(u0c  -2 a  +  6m  -  2) 

-2a  — 7 

Ul  = - 6 - 

*  COEFFICIENT  OF  -  IS  -  ^.c.2..-.?»°a2.+  2if_°a-,±  ~  2“° 

»  6 

Uo(c2  —  2a2  +  2a  —  2) 

«2  -  3g 

*  COEFFICIENT  OF  1  IS 

2uo  c3  —  3uoa2c  +  3uoac  —  3tioC  —  2a3  +  3a2  +  3a  +  108u3  —  2 

27 

2U0C3  —  3uo  arc  +  3«o<Jc  —  3uoC  —  2a3  -f  3a2  -f  3a  —  2 
U3_  _ 

*  COEFFICIENT  OF  g  IS 

c(2uoc3  —  3uoa2c  +  3uoac  —  3uoC  —  2a3  +  3a2  +  3a  —  2) 

— 

«4  IS  ARBITRARY  ! 

COMPATIBILITY  CONDITION: 

c(2uoc3  —  3uo  a2c  +  3uoac  —  3uoc  —  2a3  +  3a2  -f  3a  —  2) 

27  =  ' 

*  *  *  CONDITION  IS  NOT  SATISFIED.  *  *  * 

***  CHECK  FOR  FREE  PARAMETERS  OR  PRESENCE  OF  u0  ♦** 
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